Mapping topological order in coordinate space 
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Abstract 

The organization of the electrons in the ground state is classified by means of topological in- 
variants, defined as global properties of the wavefunction. Here we address the Chern number of 
a two-dimensional insulator and we show that the corresponding topological order can be mapped 
by means of a "topological marker", defined in r-space, and which may vary in different regions 
of the same sample. Notably, this applies equally well to periodic and open boundary conditions. 
Simulations over a model Hamiltonian validate our theory. 

PACS numbers: 73.43. Cd, 03.65.Vf, 11.30.Rd 



1 



Topological insulators are sharply distinguished from normal ones by the manner in which 
the electronic ground state is topologically "twisted" or "knotted" in k-space.^'^ But topo- 
logical order must reflect a peculiar organization of the electrons even when the concept of 
k-space does not apply, such as for inhomogeneous systems, as well as for finite systems 
within open boundary conditions. We address here the archetypical topological invariant, 
namely the first Chern number C, defined for a many-electron system in two dimensions 
(2d), and we show that the corresponding topological order also bears a very clear signature 
in r-space. We introduce a "topological marker" , which may vary in different regions of the 
same sample, and we validate our expression by means of simulations on a model Hamilto- 
nian, performed on finite samples within open boundary conditions. Our test cases include 
crystalline as well as disordered samples, and heteroj unctions. 

For a lattice-periodical system of independent electrons the Chern number (a.k.a. TKNN 
invariant^) C is expressed as a 2d Brillouin-zone integral. For a disordered and macroscopi- 
cally homogeneous system C has a known expression in a supercell framework,^'^ also formu- 
lated in k-space. The concept of k-space is rooted in the periodic boundary conditions (or 
generalizations thereof), while instead our topological marker samples the electronic ground 
state locally. The choice of boundary conditions becomes irrelevant in the limit of a large 
sample. 

For a system of independent electrons, within either periodic or open boundary conditions, 
the ground state is uniquely determined by the one-particle density matrix, a.k.a. ground- 
state projector P(r, r'); it is a "nearsighted"^"^ operator, exponentially decreasing with |r— r'| 
in insulators even when C 7^ 0.^ Our major result is expressing the topological marker in 
terms of P directly, Eqs. (8) and (9) below, where the one-particle orbitals do not appear. 

Let Un'kiy) — e~*''''"'0„k(i") be the periodic part of the Bloch orbitals, normahzed in the 
unit cell of area Af.. The standard expression of the Chern invariant in a 2d lattice-periodical 
insulator is^° 



we assume single occupancy (a.k.a. "spinless electrons") throughout. In Eq. (1) A^c is the 
number of occupied bands, and the integral is over the Brillouin zone; C is guaranteed to 
be an integer and is gauge invariant, i.e. invariant either by unitary transformations of the 
occupied orbitals among themselves, or by a change of the magnetic gauge. It characterizes 
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therefore the many-electron ground state, not the Hamiltonian. 
We start inserting a complete set of states into Eq. (1) 
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where the missing terms are real. Then, by some manipulations which are standard in 
linear-response theory, we have 

(Mn'k|VkMnk) = -^(^n'k| r |t/'nk), n 7^ n . (3) 

We stress that, while the position operator r is ill-defined within periodic boundary 

conditions, its off-diagonal elements over the Hamiltonian eigenstates are well defined; 
more accurately, Eq. (3) should be interpreted as a definition of such elements. Then 
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/o N2^"^ Xl X! /"^^ /c?k'(V'nkk|V'n'k')(V'n'k'|y|V'nk), (4) 

where the second hue owes to the fact that the matrix elements vanish for k 7^ k'. Next we 
recognize the ground-state projector P and its complement Q — 1 — P 

<5 = 7T7^ Z_\ /f^k'|V'„'k')(^n'k'|- (5) 

Eq. (4) becomes then the trace over a crystal cell of a real-space operator: 

C--i-^Im trceii{PxQy} 

TT Ac 

4:TT 

= — Im tT,,n{PxPy}, (6) 

where the symmetry of the missing term yields the second line. By exploiting the general 
properties of projectors and of the trace, Eq. (6) can be recast in several equivalent ways. 
For lattice models, a similar real-space formula has been demonstrated in 2006 by Kitaev;^^ 
to the best of our knowledge, we are providing the first proof which does not rely on lattice 
models and generalizes Kitaev's result to realistic implementations. 



Subsequent work adopting Kitaev's formula was invariably rooted in k-space within a 
toroidal geometry, for a system without boundaries, and was based on traces. ^^"^^ Finite 
systems within open boundary conditions look problematic. In fact, if we replace the trace 
over the cell with the trace over the whole sample, the identity 

Im tr{PxPy} = ^tr{ [PxP, PyP] } (7) 

guarantees a zero result, whenever P projects over a finite-dimensional manifold. This 
confirms that the global topology is trivial within open boundary conditions, and also hints 
that traces must be avoided when addressing finite and/or inhomogeneous samples. 

At variance with previous work based on Kitaev's formula, wc propose here to directly 
address the commutator in Eq. (7) before taking the trace. Let X be the projected x- 
coordinate 

X{r, r') = j dr" P{r, r")x"P{r", r'), (8) 

and similarly Y; we then identify the topological marker with the local Chern number as^'^ 

€{r) = -2ni J dr' [X{r, r')Y{r', r) - Y{r, r')X{r', r) ]. (9) 

Our definition holds within both periodic and open boundary conditions; given the short- 
sightedness of P, in a region of crystalline periodicity the cell average of €{r) coincides with 
the Chern number C owing to Eq. (6). We expect the dimensionless function €{r) to fluctu- 
ate over microscopic dimensions; in the nonperiodic case, the cell average has to be replaced 
with the macroscopic average, defined as in electrostatics (see e.g. Jackson^^). 

The gauge invariance of C(r) as defined in Eq. (9) deserves a comment. The ground- 
state projector P is invariant by unitary transformations of the occupied orbitals among 
themselves, but not by a change of the magnetic gauge. However, the unitary operator 
which transforms P is local in coordinate space, thus ensuring gauge invariance of C{r). 

We validate our formal findings by performing simulations on the Haldane model 
Hamiltonian;^^ it is comprised of a 2d honeycomb lattice with two tight-binding sites per 
primitive cell with site energies ±A, real first-neighbor hoppings ti, and complex second- 
neighbor hoppings t2e'^^'^. As a function of the parameters, this 2d model system may 
have either C = OorC = ±l, according to the phase diagram shown in Fig. 1. This 
model has been previously used in several simulations, providing invaluable insight into 
orbital magnetization^'^^'^^ as well as into nontrivial topological features of the electronic 
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FIG. 1. Chern number of the bottom band of the Haldane model as a function of the parameters ^ 
and A/t2 {h = l,t2 = 1/3). The points marked with letters (a-e) in this phase diagram are relevant 
for the subsequent discussion and figures. In order to avoid special features the parameter is not 
a multiple of 7r/4. 

wavefunction.^'^'^^'^^'^^ At half filling the system is insulating, except when A — t2sm(f) — 0. 
In the present work we study, within open boundary conditions, finite flakes of rectangular 
shape cut from the bulk, as shown in Fig. 2. Wc have addressed homogenous samples 
where the Hamiltonian is chosen from various points of the phase diagram. Fig. 1, as well 
as disordered and inhomogeneous samples. 

Two typical plots for crystalline samples arc shown in Fig. 3, where we have chosen the 
two points (b) and (c) in Fig. 1, with C = and C = 1, respectively. The plots confirm 
that the local Chern numbers €{i) are equal to either or 1 (as expected) in the bulk of 
the sample, while they deviate in the boundary region. In both cases the negative values 
compensate for the positive ones, given that the sum of the over the whole sample 
vanishes. This compensation is most interesting when C — 1 (right panel). A size analysis 
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FIG. 2. A typical flake, with 2550 sites, showing the honeycomb lattice of the Haldane model^^. 
The 50 sites on the horizontal line will be used in all the subsequent one-dimensional plots. Black 
and grey circles indicate nonequivalent sites (with onsite energies ibA) 



shows that the minimum negative value scales like L (linear dimension of the sample): 
the reason is that the number of bulk sites scales as L^, while the perimeter scales as L. 

We have studied both polar (A ^ 0) and nonpolar (A = 0) cases. While in the latter case 
the two sites are equivalent, they are no longer so in the former case. This clearly appears 
in the site occupancies, also shown Fig. 3. What is surprising, is that the corresponding 
€{i) values do not show any site alternance, while we expect only their cell (or macroscopic) 

average to be equal to one. We conjecture this to be due to extra symmetry present in the 
Haldane model Hamiltonian, actually broken in disordered samples, discussed below (see 
Fig. 5). 

We have also investigated a few points in the phase diagram close to the transition between 
C — and C = 1 at fixed A/t2 — 3.67 and various cj) values. Given the finite size of the 
system the transition cannot be sharp. The exact transition for an infinite system occurs at 
(p/n = 0.25; our results show that in the bulk of the sample the local Chern number is zero 
up to 0/7r ~ 0.17 and one from (p/ir ~ 0.29 onwards. At intermediate values the boundary 
region broadens considerably and indeed invades the whole sample: this is shown in Fig. 4. 

Typical results for disordered — and macroscopically homogenous — samples are shown in 
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FIG. 3. Local Chern number (top) and site occupancy (bottom), for the 50 sites along the line 
shown in Fig. 2. Left panel: point (b) in the phase diagram, Fig. L Right panel: point (c). Notice 
the different scales. 



Fig. 5. In the left panel the sign of A alternates between the two sublattices, while its 
modulus is chosen at random (with uniform distribution) in the (a-b) segment of Fig. 1. In 
the right panel the value of A is chosen at random in the (c-d) segment. It appears clearly 
that the local Chern numbers in the bulk of the sample oscillate around a macroscopic 
average C — (left panel) and C — 1 (right panel). 

Next we show in Fig. 6 our topological marker across an hetero junction between regions 
of different topological order, in two typical cases: a normal insulator joined to a C = 1 
insulator, and a junction where C changes sign. In both cases the marker maps very per- 
spicuously the actual topological order in the two bulklike regions, while it oscillates at the 
interface and at the sample boundary. The virtue of our r-space approach is clearly demon- 
strated; the conventional k-space approach to topological order cannot separate different 
regions of inhomogeneous samples. 

Finally we analyze the present results from the viewpoint of the modern theory of the 
insulating state. ^'^ Both Eqs. (1) and (6) look hke the imaginary offdiagonal part of a more 
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FIG. 4. Local Chern number for a few points on the line A/f2 = 3.67, i.e. on the (b-c) segment in 
Fig. 1, close to the transition from C = to C = 1. The exact transition occurs at ^/tt = 0.25; 
our five plots correspond (bottom to top) to ^jv: = 0.17,0.25,0.27,0.29,0.33. 



general tensor; its corresponding symmetric part is real and measures indeed the localization 
of the electronic ground state in any homogeneous insulator. The key ingredient is the 
localization tensor (rQr^)c, a.k.a. second cumulant moment of the electron distribution 
(Greek subscripts are Cartesian indices); it has the dimensions of a squared length and its 
trace provides the gauge-invariant part of the quadratic spread of the Wannier functions, 
according to the Marzari-Vanderbilt theory. Notice that localized Wannier functions do 
not exist whenever C ^ Of^ nonetheless {rarp)c remains well defined and finite in any 
insulator. ^'^^ 

The direct link between {rarp)^ and C has been investigated elsewhere within periodic 
boundary conditions;*'^^ to see the relationship to Eq. (6) we write the localization tensor 
as7'27 

{rarp)c = jj-t^ceii{raPrpQ}, (10) 

whence C = 47rnoIm (rir2)c {no — Nc/A^ is the density). It has been shown''"'^ that Eq. (10) 
generalizes to finite systems within open boundary conditions, just taking the trace over the 
whole system and dividing by the total number of electrons. The tradeoff is that (rQ,r^)c 
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FIG. 5. Local Chern number (top) and site occupancy (bottom), for disordered systems (see text). 
Left panel: disordered system along the line (a-b) in the phase diagram, Fig. 1. Right panel: line 
(c-d). Notice the different scales. 



becomes then real symmetric, in agreement with the present findings. 

In conclusion, we have addressed a system of independent spinless electrons in 2d, whose 
topological order is classified by means of the archetypical topological invariant: the Chern 
number C. We have found the explicit form of a local Chern number £(r). It is a gauge- 
invariant microscopic function, whose macroscopic average coincides with C in crystalline 
samples. Notably, the boundary conditions (either periodic or open) are irrelevant in the 
definition of €{r). For disordered and/or inhomogeneous samples the macroscopic average 
of €{r) is a marker which detects the kind of topological order in any macroscopically 
homogeneous region: for instance either in a disordered sample, or across an hetero junction. 

At the root of our local description of topological order is the "nearsightedness" of the 
ground-state density matrix. Since this is a very general feature of insulators, it is possi- 
ble that any kind of topological order^'^ — described by invariants other than C — could be 
addressed via the appropriate local marker. 
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FIG. 6. Local Chern number (top) and site occupancy (bottom), across hetero junctions. Left 
panel: Hamiltonian parameters as in (a) and in (b) for left and the right halves of the sample, 
respectively. Right panel: parameters as in (e) and in (c) for left and the right halves of the sample. 
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